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CN ■ Abstract 

Q_i! We study the direct CP violation in B° — > p°(uj)p (uj) — > tt + tt~ tt + tt~ (with unpolarized p°(uj)) 

< 

via the p — lj mixing mechanism which causes a large strong phase difference and consequently 
a large CP violating asymmetry when the masses of the ir + ir~ pairs are in the vicinity of the u 
resonance. Since there are two p(ui) mesons in the intermediate state p — to mixing contributes 

On' 

twice to the first order of isospin violation, leading to an even larger CP violating asymmetry 
(could be 30% - 50% larger) than in the case where only one p(uj) meson is involved. The CP 
violating asymmetry depends on the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements and 
the hadronic matrix elements. The factorization approach is applied in the calculation of the 



o 



X 



hadronic matrix elements with the nonfactorizable effects being included effectively in an effective 
parameter, A^ c . We give the constraint on the range of -/V c from the latest experimental data for the 



CO 

branching ratios for B° — > p°p° and B° — > p + p . We find that the CP violating asymmetry could 



be very large (even more than 90% for some values of N c ). It is shown that the sensitivity of the 
CP violating asymmetry to iV c is large compared with its smaller sensitivity to the CKM matrix 
elements. We also discuss the possibility to remove the mod (tt) ambiguity in the determination of 
the CP violating phase angle a through the measurement of the CP violating asymmetry in the 
decay B° — > p° (to) p° (u) — > ir + ir~ tt + ir~ . 
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I. INTRODUCTION 



Although CP violation has been a central concern in particle physics since it was first 
observed in the neutral kaon system more than four decades ago [1] the dynamical origin of 
CP violation still remains an open problem. CP violation in the framework of the Standard 
Model (SM) is supposed to arise from a weak complex phase in the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix which is based on quark flavor mixing 2, jij]. Therefore, the study 
of CP violation is essential to the test of the CKM mechanism in the SM. 

Besides the kaon system much more studies have been carried out on CP violation in 
the B meson system both theoretically and experimentally in the past few years. It was 
suggested theoretically that large CP violating asymmetries should be observed in the ex- 
periments for B mesons [4j. This important prediction has already been confirmed by the 
experiments of BaBar and Belle etc. through the measurements on CP violation in several 

n 

decay channels of B mesons such as B° -> J/ipK° s and B° -> K + tt~ |5(. From the summer of 
2007, the Large Hadron Collider (LHC) at CERN will start to contribute to the exploration 
of CP violation in the B meson system in a more accurate way due to its much higher 
statistics. This will also provide an opportunity to discover new physics beyond the SM. 

In the decay process we have the so-called direct CP violation which occurs through 
the interference of two amplitudes with different weak phases and strong phases. The weak 
phase difference is directly determined by the CKM matrix. On the contrary, the strong 
phase is usually due to complicated strong interaction and hence difficult to control. Since 
a large strong phase difference is required for a large CP asymmetry, one needs to appeal to 
some phenomeno logical mechanism to get such a large strong phase difference. The charge 
asymmetry violating mixing between p° and uj (p — uj mixing) has been applied for this 
purpose in the past few years. From a series of studies for CP violation in some decay 
channels of heavy hadrons including B, and D, it has been found that p — uj mixing can 
provide a very large strong phase difference (usually 90 degrees) when the mass of the decay 
product of p(uj), 7r + 7r~, is in the vicinity of the uj resonance Q, si, 1(3] . Furthermore, it 
has been shown that the measurement of the CP violating asymmetry for these decays can 
be used to remove the mod (tt) ambiguity in the determination of the CP violating phase 
angle a. 

In this paper, we will investigate the CP violating asymmetry for the decay B° — > 
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p°(uj)p°(uj) — > 7r + 7r _ 7r + 7r _ . This process is unique since it has two p(uj) mesons in the 
intermediate state, each of them contributing p — uj mixing. One can expect that there 
should be a bigger CP violating asymmetry than in the case where p — uj mixing only 
contributes once. It will be shown from our explicit calculations that this is true indeed. 
The CP violating asymmetry in the case of double p — uj mixing could be 30 - 50% bigger 
than that in the case of single p — uj mixing, depending on the value of N c and q 2 jm\ (see 
the meaning of N c and q 2 jm\ below). 

In our calculations of the CP violating asymmetry, hadronic matrix elements for both 
tree and penguin operators in the effective Hamiltonian are involved. These matrix elements 
are controlled by the effects of nonperturbative QCD which are difficult to handle. In order 
to extract the strong phase difference we will use the factorization approximation, in which 
one of the currents in the Hamiltonian is factorized out and generates a meson, assuming 
the vacuum intermediate state saturation. In this way, the decay amplitude becomes the 
product of two matrix elements. Such factorization scheme was first argued to be plausible 
in energetic decays like bottom-hadron decays 111] 121], then was proved to be the leading 
order result in the framework of QCD factorization when the radiative QCD corrections of 
order 0{a s {m^)) (m& is the b-quark mass) and the 0(l/m&) corrections in the heavy quark 
effective theory are neglected 13j . Since the nonfactorizable contributions are ignored in the 
factorization scheme we introduce an effective parameter, N c , in order to take into account 
nonfactorizable contributions effectively. In this way, the value of N c is not the color number 
(3) any more, but should be determined by experimental data. In the present work, this 
will be done by comparing the theoretical results with the experimental data for the decay 
branching ratios for the processes B° — > p°p° and B° — > p + p~ . 

The remainder of this paper is organized as follows. In Sec. [IT] we briefly present the 
effective Hamiltonian, the Wilson coefficients and the CKM matrix elments. In Sec. IIHI we 
give the formalism for the CP violating asymmetry in B° — > p°(uj)p°(uj) — > 7t + 7t~7t + tt~ via 
p — uj mixing. Then we give the calculation details of the strong phase difference and the 
numerical results for the CP violating asymmetry. In Sec. IIV1 we calculate the branching 
ratios for B° — > p°p° and B° — > p + p~ and present the range of N c allowed by the latest 
experimental data for these decays. In the last section, we give a summary and discussion. 
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II. THE EFFECTIVE HAMILTONIAN AND THE CKM MATRIX 



In order to calculate the direct CP violating asymmetry one needs to use the following 
effective weak Hamiltonian based on the operator product expansion [l^ : 

Q 10 

H± B =i = [ E VrtKfaoz + c 2 o u 2 ) - v tb v; q £ ctOt] + h.c, (i) 

q=d,s i=3 

where q(i=l,...,10) are the Wilson coefficients, V u b, V uq , V t b and V tq are the CKM matrix 
elements. The operators have the following form: 

Oi = q a 1^(l - 7 5 )w / 3M / 37 / "(l - 75)&a, 
0^ = ^(1-75)^(1-75)6, 
3 = ?7m(1 - 75)6^g'7 M (l - 7s)g', 

q' 

Oa = q a J^(l - 75)6/3 ^7 M (1 - 75)g a , 

5 = ?7m(1 - 75)6^g'7 M (l + 7 5 )g, 

6 = g« TM (l - 75)6/3 E ^7 M (1 + ^Va, 

7 = 2^(1 - 75)6 Zq'V'^i 1 + TS)?', 

8 = 2^7m(1 - 7s)6/3 22 e 9'^/37 M (l + 75)^, 

9' 

09 = 2^(1 - 75)6^e g /g'7^(l - 7 5 )g', 

q' 

O w = 2^7^(1 - 75)6/3 2^ <V0 VC 1 - 7sK, (2) 

<?' 

where a and /3 are color indices, and q' = u, d or s quarks. In Eq. (T5]) O" and 0% are tree 
operators, 3 -0 6 are QCD penguin operators, and 7 -O w arise from electroweak penguin 
diagrams. 

The Wilson coefficients are known to the next-to-leading logarithmic order [14]. fis|. They 
are renormalization scheme dependent since the renormalization prescription involves an 
arbitrariness in the finite parts in the renormalization procedure. The physical quantities 
should be renormalization scheme independent. Since the radiative QCD corrections are 
not included in the factorization approach we work, the hadronic matrix elements do not 



carry any information about the renormalization scheme dependence 1 . Therefore, we choose 
to use the renormalization scheme independent Wilson coefficients which are defined in 
Refs. Q, 
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17l | so that the CP violating asymmetry we obtain is renormalization scheme 
independent. The renormalization scale \i is chosen as the energy scale in the decays of 
the B meson, 0{m b ). When fi = 5 GeV, these renormalization scheme independent Wilson 



coefficients take the following values [16j, 117] : 

ci = -0.3125, c 2 = 1.1502, 
c 3 = 0.0174, c 4 = -0.0373, 
c 5 = 0.0104, c 6 = -0.0459, 
c 7 = -1.050 x 10~ 5 , c 8 = 3.839 x 10~ 4 , 

c 9 = -0.0101, cio = 1-959 x 10~ 3 . (3) 

The matrix elements of the operators should be renormalized to the one-loop order. 
This results in the effective Wilson coefficients, 4, which satisfy the constraint 

c i (m b )(O r (m b )) = c' l (O l ) t * ee , (4) 

where (0,) tree are the matrix elements at the tree level, which will be evaluated in the 

nn 

factorization approach. From Eq. (jlj), the relations between c[ and Cj are |16l. Il7| 

c 'l = c l; c 2 = C 2) 

c 3 = c 3 — P s /3, c 4 = C4 + Pg, 
c' 5 = c 5 - P s /3, c' 6 = c 6 + P s , 

Cy = C7 + P e , Cg = Cg, 

4 = C 9 + P e , 4o = c 10, (5) 

where 

P s = (« s /8vr)c 2 [10/9 + G(m c , fi, q% 

P e = (a em /97r)(3 Cl + c 2 )[10/9 + G(m c ,/i,g 2 )], 



1 It has been shown that in the QCD factoization approach the renormalization scheme dependence of the 
Wilson coefficients and that of the hadronic matrix elements cancel 13( . 
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G(m c , fj,, q 2 ) = 4 / dxx(x — l)ln 



— x(l — x)q 2 



with 



'o A* 
where m c is the c-quark mass and g 2 is the typical momentum transfer of the gluon or photon 
in the penguin diagrams. G(m c , /i,q 2 ) has the following explicit expression 181 ]: 



3mG 



i 

3 \ fi 2 3 q 2 \ q 2 



2 I . mi 5 ,m; / mt 
- In— ^ - - - + 1 + 2-£ 
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^ l-Jl-4^ 



r( i+ w- 4 v 



(6) 



The value of g 2 is chosen to be in the range 0.3 < q 2 ' jm\ < 0.5 . From Eqs. (j3J) (J5]) 
(|SJ) we can obtain numerical values of d i which are listed in Table [U where we have taken 
cv s (7«z)=0.118, a em (m b )=l/132.2, m b =5 GeV, and m c =1.35 GeV. 



TABLE I: Effective Wilson coefficients for the tree operators, electroweak and QCD penguin op- 
erators 0, Q 



< 


g 2 /mg=0.3 




q 2 /ml=0.5 




-0.3125 




-0.3125 




1.1502 




1.1502 


4 


2.433 x 10" 2 + 1.543 x 10" 


H 


2.120 x 10" 2 + 5.174 x 10~ 3 f 




-5.808 x 10" 2 - 4.628 x 10 


-h 


-4.869 x 10~ 2 - 1.552 x 10^ 2 i 




1.733 x 10~ 2 + 1.543 x 10 _ 


H 


1.420 x 10~ 2 + 5.174 x 10 3 i 




-6.668 x 10~ 2 - 4.628 x 10 


~ 3 i 


-5.729 x 10~ 2 - 1.552 x 10^ 2 i 




-1.435 x 10~ 4 - 2.963 x 10 


-H 


-8.340 x 10" 5 - 9.938 x 10~ 5 i 




3.839 x 10~ 4 




3.839 x 10~ 4 




-1.023 x 10~ 2 - 2.963 x 10 


-H 


-1.017 x 10~ 2 - 9.938 x 10~ 5 i 


c 10 


1.959 x 10~ 3 




1.959 x 10~ 3 



The CKM matrix, which should be determined from experiments, can be expressed in 
terms of the Wolfenstein parameters, A, A, p and rj 191 ]: 
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/ 1 



|A 2 
-A 



A AX 3 (p - irj) 
1 - \\ 2 AX 2 



\ 



(7) 



yA\ 3 (l- p-irj) -AX 2 1 j 

where 0(A 4 ) corrections are neglected. The latest values for the parameters in the CKM 



matrix are 



20j: 



A = 0.2272 ± 0.0010, A = 0.818^;^, 
n - 2 21 + 064 ft - S40+ - 017 



where 



A 2 A 2 
P = p(l-y)> *7 = *7(1 - y)- 



From Eqs. ([HD ( [9} we have 



0.198 < p < 0.293, 0.302 < rj < 0.366. 



(8) 



(9) 



(10) 



III. CP VIOLATION IN B° -» p°(u)p°(uj) -» vr+vr-vr+vr- 



Formalism 



Letting A (A) be the amplitude for the decay B° — > 7r + 7r 7r + 7r (S° — > 7r + 7r 7r + 7r ) one 



has: 



A : 



Tr+Tr-Tr+Tr^l^lS ) + (7r + 7r"~7r + 7r \H P \B°), 
yn-7i + 7i-\H T \B°) + (n + n-n + n-\H p \B ), 



(11) 
(12) 



with if T and if p being the Hamiltonian for the tree and penguin operators, respectively. 

We can define the relative magnitude and phases between the tree and penguin operator 
contributions as follows: 



A = (n + n-n+n-\H T \B°)[l + T S s +*>\ 
A = (n + n-n + n-\H T \B°)[l + re l{5 '^ 



(13) 
(14) 



where 6 and <fi are strong and weak relative phases, respectively. The phase <fi can be 
expressed as a combination of the CKM matrix elements: 4> = ^g[(V t bV t * d ) / (V u bV* d )]. As a 



result, sin0 is equal to sina with a being defined in the standard way 20j. The parameter 
r is the absolute value of the ratio of penguin and tree amplitudes: 



(15) 



(tT + 1T 7T + 7T 


H p 


B°) 




H T \B°) 



The CP violating asymmetry, a, can be written as 

— 2rsin5sin0 



\A 


2 - 


A\ 2 


\A 


2 + 


A 


2 



(16) 



1 + 2rcos^cos0 + r 2 

In order to obtain a large signal for direct CP violation, we need some mechanism to make 
sin<5 large. It has been found that p — uj mixing has the dual advantages that it leads to a 
large strong phase difference and is well known , [h]]. With this mechanism, to the 

first order of isospin violation, we have the following results when the invariant masses of 
7r + 7r~ pairs are near the uj resonance mass: 



Tr+TT-TT+Tr-I^IE ) 
H P \B°) 



H ~ gl 
p n t + p 



IT IT 11 IT 



2 PPi 
P 
2 



sis, 



"^-pujPpuj nPpP' 



(17) 
(18) 



Here t pp (p pp ) and t puJ (p pu) ) are the tree (penguin) amplitudes for B — > p°p° and B° — > 
respectively, g p is the coupling for p° — > tt + tt~, T\. puJ is the effective p — u mixing amplitude 
which also effectively includes the direct coupling uj — > tt + tt~, and sv(V=p or uj) is the 
inverse propagator of the vector meson V, 



s v = s - my + im v T Vl 



(19) 



with y/s being the invariant masses of the tt + tt pairs (we let the invariant masses of the 



two tt + tt pairs be the same). Eqs. 

fin 

single p — uj mixing is involved [8|, [9|, 



T71) (fT8l) have different forms from the case where only 
101 ]: there is a factor of 2 in front of the effective p — uj 



mixing amplitude, n pa; , since p — uj mixing contributes twice to the first order of isospin 
violation. Furthermore, we have g p and s 2 p instead of g p and s p as before due to two p — > tttt 
couplings and two p propagators (note that s 2 term is of the second order of isospin violation 
and hence is ignored). 

As mentioned before, the direct coupling uj — > it + it~ has been effectively absorbed into 
U puJ [21]. This leads to the explicit s dependence of H puJ . In practice, however, the s 
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dependence of n pw is negligible. Making the expansion U piJJ (s) = IL^fm, 2 , ) + (■§— m„,)U ( m 2 ), 
the p — uj mixing parameters were determined in the fit of Gardner and O'Connell 221 ] : 

y\zU pu (ml) = -3500 ± 300MeV 2 , 
Ml^m 2 ) = -300 ± 300MeV 2 , 

n^(m 2 ) = 0.03 ±0.04. (20) 

From Eqs. (jTTj) (fl3l (jTTj) ([15]) one has 



2hlp a ;tp a ; -|- S^tpp 

where the factor of 2 in front of U puJ arises from the involvement of double p — uj mixing. 
Defining 

Pp« = r / e i&+0 V= ae ^ Ppp.^^^ (22) 
^pp ^pp Ppii) 

where <5 a , 5p and <5 g are strong phases, one finds the following expression from Eqs. (T2"Tj) f[2"2"l) : 

re l5 = r'e^^^-^. (23) 
In order to obtain the CP violating asymmetry in Eq. ffTB]) . sin0 and cos</> are needed, where 



191 ] . one 



<j) is determined by the CKM matrix elements. In the Wolfenstein parametrization 
has 

sin = 2 (24) 

cos0= . (25) 

V[p(!-P) -?? 2 ] 2 + r, 2 

B. Calculational details 

With the Hamiltonian given in Eq. ([1]) we can evaluate the matrix elements for B° — > 
p° (uj) p° (uj) . In the factorization approximation, p°(uj) is generated by one current which has 
the appropriate quantum numbers in the Hamiltonian. For this decay process, the amplitude 
can be written as the product of two matrix elements after factorization, i.e. (omitting 
Dirac matrices and color labels): (p°(uj)\(qq)\0)(p°(uj)\(db)\B ) (q = u,d), where (qq) and 
(db) denote the V — A currents, <?7p(l — 7s)g and dj^l — 75)6, respectively. Since p° and 



u are vector mesons the amplitude for B° — > p°(uj)p°(uj) may be polarized or unpolarized. 
Here we investigate the later case. Defining 

(p°\(uu)\0)(p \(db)\B )=T, (26) 

one has 

T =-(p°\(dd)\0)(p \(db)\B ) 
= -(p \(u U )\0}(u;\(db)\B } 
= (p°\(dd)\0)(u J \(db)\B ) 

= (u\(uu)\0}(p°\(db)\B }. (27) 
After factorization, the contribution to t pp from the tree level operator 0\ is 

(p°p |O^|fi°) = 2(p°\(uu)\0}(p°\(db)\B } = 2T. (28) 
Using the Fierz transformation the contribution of is (1/N C )T. Hence we have 

t pp = 2(c[ + ^-AT. (29) 



It should be noted that in Eq. (T291 we have neglected the color-octet contribution which 
is nonfactorizable and difficult to calculate. Therefore, N c should be treated as an effective 
parameter and may deviate from the naive value 3 . In the same way we find that 

tpuj = 0. This lead to 

ae iSa = 0, (30) 



from Eq. 

In a similar way, we can evaluate the penguin operator contributions p pp and p puJ with 
the aid of the Fierz identities. From Eq. (T221 we have 



Ufi _ -2(<4 + wA) +3(4 + ^j) + 3(4 + ^c' 10 ) + ( c ; + ^c; 

2 (4 + ^4) + 2 + ^3) + 2 (4 + ^c' 6 ) - (4 + ^c' 8 ) - (c' 9 + i4 ) - (c' 10 + ^c' 9 

(31) 

-2(4 + ±d) -2(d + 4r&) -2(^4 + ^4^1 + fc' + 4r<&) + (d + ttC>) + (c' + ±-d 



2 ( c 'i + ^ c 4 



x £j& , (32) 
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where 



v ub v: d 



(33) 



Numerical results 



We have several parameters in the numerical calculations: q 2 , N c , and the CKM matrix 
elements. Since N c includes nonfactorizable effects, which cannot be evaluated accurately at 
present, we choose to treat it as a parameter and determine its range from the experimental 
data. Then, we can extract an allowed range for N c from a comparison of the theoretical 
results and the experimental data. By doing this, we get the range of N c as 2.74(2.81) < 
N c < 4.77(4.92) for q 2 jm\ = 0.3(0.5). This will be discussed in detail in Sec. HYl The most 
uncertainties due to the CKM matrix elements come from p and r\ since A is well determined 
(see Eq. (JS}) and since the CP violating asymmetry is independent of the Wolfenstein 
parameter A. Therefore, in our numerical calculations, we take the central value for A and 
only let (p, 77) vary between the limiting values (p min , r] min ) and (p max , r) max ). In fact, explicit 
numerical results show that the CP violating asymmetry is very insensitive to A. 

In the numerical calculations, it is found that for a fixed N c there is a maximum value, 
dmaxi for the CP violating parameter, a, when the invariant masses of the 7r + 7r~ pairs are 
in the vicinity of the uj resonance. This is shown explicitly in Fig. 1. For q 2 /m 2 = 0.3(0.5) 
and N c = 2.74(2.81), the maximum CP violating asymmetry varies from around -91.1% 
(-70.1%) to around -96.1% (-77.8%) as (p,rj) change from (p max ,r] max ) to (p m in,Vmin); For 
q 2 /m 2 = 0.3(0.5) and N c = 4.77(4.92), the maximum CP violating asymmetry varies from 
around 55.8% (28.9%) to around 53.2% (22.9%) when (p, rj) change from (p ma x,Vmax) to 

sprain j T]min) ■ 

Our results show that the p — uj mixing mechanism produces a large sin<5 in the allowed 
range of N c , which is necessary for a large CP violating asymmetry. The involvement of 
double p — uj mixing in B° — > p°(uj)p°(uj) — > 7r + 7r~7r + 7r~ gives rise to a factor of 2 in Eq. 
([211 in front of H puJ . This makes the CP asymmetry even larger than in the case where the 
p — uj mixing contributes only once. Fig. 2 shows explicitly the comparison between these 
two cases for iV c = 2.74(2.81) and q 2 /m 2 = 0.3(0.5). The maximum asymmetry with the 
involvement of single p — uj mixing, for q 2 /m 2 = 0.3(0.5) and N c = 2.74(2.81), is around 
-55.2% (-20.0%) for the set (p ma x,Vmax) and -63.2% (-23.8%) for the set {p min ri min ). For 
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q^/ml = 0.3(0.5) and iV c = 4.77(4.92), we find that a max is around 26.5% (-1.49%) for the 
set (pmax,Vmax) and 25.1% (-1.50%) for the set (p m inVmin) in the case of single p — to mixing. 

The reason that double p — uj mixing leads to a larger CP violating asymmetry than 
in the case of single p — uj mixing is that sinS becomes bigger in the case of double p — uj 
mixing than in the case of single p — uj mixing. This can be seen explicitly from Fig. 3. The 
involvement of double p — uo mixing may also change the value of r as can be seen from Fig. 

4. However, as found from our detailed analysis for the influence of r on the CP violating 
asymmetry, the effect of the change of r on a is small compared with the change of sin^ due 
to the involvement of double p — uj mixing. 

It is noted that when N c is around 2.81 and 4.92 in the case q 2 jm\ = 0.5, we could also 
have large CP asymmetries when y/s is far away from the uj resonance for all the allowed 
values of the CKM matrix elements (see Figs. 1(b) and 5(b)). In these cases, the effective 
p — uj mixing contributes little and the large CP asymmetry is caused by the effective Wilson 
coefficients, which can also give a large strong phase, S, since they are complex numbers. 

In most direct CP violating decays such as B° — ► p°(u;)7r — ► n + n~ii [9|, llO] and some 
other processes, the involvement of p — uj mixing leads to the result that the strong phase, 

5, passes through 90° (sin<5=l) at the uj resonance. However, in the decay we are discussing 
this does not happen in the allowed range of N c . Instead, the absolute value of sin<5 just 
gets close to 1, but does not equal 1 (see Fig. 3), even though it is enough to give large CP 
asymmetry, especially when double p — uj mixing is involved. 

Figs. 5 and 6 show the dependence of the CP violating asymmetry and sin#, respectively, 
on both N c and y/s. One can see that the CP asymmetry strongly depends on iV c . Take 
Fig. 5(a) as an example (for q 2 ' jm\ = 0.3 and maximum p and 77): when N c < 3.68, one 
gets minus asymmetry around the uj resonance, whereas when N c > 3.68 the CP violating 
asymmetry becomes positive. 

It can be seen from Fig. 5 that when N c takes the critical value, — c^/c^ ~ 3.68, the CP 
violating asymmetry becomes zero. This is because t pp = at this point (as can be seen 
from Eq. ( 1291) easily) and hence the penguin operator contributions dominate. Furthermore, 
the sign of sin5 and hence the sign of the CP violating asymmetry change at this point. It 
would be interesting to see whether or not iV c can take this value in the future when more 
accurate experimental date are available. From most previous studies, it seems that N c is 



usually less than this critical value [7, 8, |9|, llC)|, l23|. If this is true for B° — > p (uj)p°(uj), then 



12 



the sign of sin<5 would remain unchanged. Then, one could remove the mod (jr) ambiguity 
in the determination of the CP violating phase angle a (through sin2a) by measuring the 
CP violating asymmetry in B° — > p°{uj)p (uj) — > 7r + 7r~7r + 7r~. 



IV. BRANCHING RATIOS FOR B° -> p°p° AND 5° p+p~ 



A. Formalism 



If the decay amplitude for B — > V{V 2 (Vi,V 2 denote vector mesons) has the form 
A(B -> ViV 2 ) = aX {BViy2 \ where X {BViy '^ denotes the factorizable amplitude with the 



form (V2 1 (^2^3) |0) (Vi I (gi6) I i3) , then the decay rate is given by [23j 



Y{B - V 1 V 2 ) = -^-\ a (m B + mi )m 2 f V2 Af v ^m 2 )\ 2 H, (34) 



tfrm 2 B 



where a is related to the CKM matrix elements and Wilson coefficients, fv 2 is the decay 
constant of V 2 , p c is the cm. momentum of the decay particles, and mi(m 2 ) are the 
masses of the B meson and the vector meson Vi(T^), respectively, and 



where 



H= (a-bx) 2 + 2(1- c 2 y 2 ), (35) 



m 2 B — m\ — m\ 2m 2 B p 2 2rriBp c 



2m,\m 2 mim 2 (mB + mi) 21 (mg + 



A BV Hml) V BV H 

X = ,BVw y 



m 



A BV ^m 2 2 y y A BV h 



m% 



_ \/Wb - ( m i + m a) 2 ] [m B - (mi - m 2 ) 2 } 
c 2m b 

A x 1 , A BVl and V BVl in Eqs. (I35I) and (|36|) are the form factors associated with B — > Vi 
transition. 

The decay amplitudes for 5° — > p°p°, -B — > p°co> and _B° — ► p + p~ are 

A (B° - p°p°) = ai X^°'P°\ (37) 
A(5° -> p°cu) = a 2 X^ Bp °' u \ (38) 

and 

-> p+p-) = a 3 X^ + ^\ (39) 
13 



where 



ax = -j=[2aiV ub V* d - (-2a 4 + 3a 7 + 3a 9 + a 10 )V tb V t * d ] 



ol% = p(2a 3 + 2a 4 + 2a 5 — a 7 



V2 



aio)V tb V t 



td- 



G 



V2 

with cii (i = 1, 2, • • • , 10) being defined as: 



tdb 



(40) 
(41) 
(42) 



a#-i = 4,-1 + tt> for 3 = h 2, 

vv r . 



(43) 



When we calculate the branching ratios we should take into account the p — u mixing 
contribution for consistency since we are working to the first order of isospin violation. Then, 
we obtain the branching ratio for B° — > p°p°: 

2 



BR(B° - p°p°" 



Pc 



87rm B T B o 



Ot\ + ot 2 - 



2n 



(m B + m p o)m p o/ p oAi(m^ 



( 



H. 



4) 



For B° — ► p + p , we have 

5i?(5° -> p+p") 



Pc 



87rm^r £0 



(a 3 (m B + m p +)m p -fp-A 1 (m + )\ H. 



(45) 



B. Form factor models 



The form factors Ai(k 2 ), A 2 (k 2 ) and V(k 2 ) depend on the inner structure of hadrons 
and consequently depend on the phenomenological models for hadronic wave functions. We 
adopt the followin g fo rm factor models: 

Model 1(2) [24], I2J: 



V{k 2 



V(0) 



-,A 1 (k 2 



Mo) 



-,A 2 (k 2 



\-k 2 j{m\_y iV ' l-k 2 /(m 2 1+ y 
where 1/(0) = 0.33(0.395), A^O) = A 2 (0) = 0.28(0.345), m x - 
5.71GeV. 

Model 3(4) [23, bJ, bd] : 

^(0) „ ,, 2 , A 1 (0) 



V(k 



[1 - k 2 /(m{_)Y 



A^k 2 



l-k 2 /(m 2 Y 



A 2 {¥ 



AM 



l-k 2 /{m\ + y 



(46) 



5.32GeV, and min 



AM 



[1 - k 2 /{m? 1+ )] 2 ' 



(47) 
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where the form factors have double pole dependence and the parameters take the same 
values as in Models 1 and 2. 
Model 5 



26j: 



for V(k 2 



for A(k 2 ) (i=l, 2): 



V(k 2 



V(0) 



Mk 



1 — k 2 /mv)[l — aik 2 /my + a 2 k A /m 

Mo) 



41' 
V\ 



(48) 



(49) 



1 — a\k 2 /m 2 } + a 2 k A /m v '' 
where m y = m B * = 5.32GeV; V(0) = 0.31, Oi = 0.59 and a 2 = for 1/(A; 2 ); A x (0) = 0.26 
(7i = 0.73 and 0-0 = 0.10 for A X (A; 2 ); and A 2 (0) = 0.24, a x = 1.40 and o 2 = 0.50 for A 2 (k 2 ). 



Model 6 
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the form factors Ai(k 2 ), A 2 (k 2 ) andV(k 2 ) have the same form: 

/(0) 



1 — a F k 2 /m 2 B + b F k 4 /m B 



4 ' 



(50) 



where / could be At, A2, or V. The parameters /(0), and 6^ for various form factors 
are: for A u A^O) = 0.261, a F = 0.29, b F = -0.415; for A 2 , A 2 (0) = 0.223, a F = 0.93, 
6 F = -0.092; and for V, 1/(0) = 0.338, a F = 1.37, b F = 0.315. 



C. Numerical results 



As mentioned before, N c includes the nonfactorizable effects effectively, which cannot be 
handled well at present. Therefore, we treat iV c as a parameter to be determined by exper- 
imental data. Usually iV c is assumed to be universal for all decay channels in the factoriza- 
tion approach. However, it certainly could be different for different channels. Therefore, we 
choose to determine the range of iV c for B° — > p°(u)p°(uj) from the experimental data for the 
branching ratios for the decays B° — > p°p° and B° — ► p + p _ (we expect the nonfactorizable 
contributions and hence the values of iV c in these two channels are the same if isospin vio- 
lation is ignored). In order to find the range allowed for iV c we use the latest experimental 
data for the branching ratios for the two decay channels, B° — > p°p° and B° — > p + p~ [20I ]: 

BR(B° -> p°p°) < 1.1 x 10~ 6 , 
BR(i?° — > p + p~) = (2.5 ± 0.4) x 10~ 5 . (51) 
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We calculate the branching ratios for B° — > p°p° and 5° — > p + p~ with the formulae given 
in Eqs. (1341) - fj45|) in all the models for the weak form factors associated with B° — > p° 
and B° — > p + (p~), which are mentioned in the previous subsection. In addition to the 
dependence on N c , these two branching ratios also depend on the CKM matrix elements 
which are parameterized by A, A, p and rj, with the experimental values of them being given 
in Eqs. (JSJ) (TTOl) . Since each of these parameters has some uncertainty, we let each of them 
vary in its allowed range when we calculate the branching ratios. Then, for each set of the 
values for the parameters A, A, p and 77, we obtain a range of N c which is allowed by the 
experimental data for the branching ratios for both B° — > p°p° and B° — > p + p~. This is 
shown in Figs. 7 and 8 for a special set of CKM matrix parameters when q 2 /m 2 = 0.3. 
Repeating this process for various sets of the values for A, A, p and r\ and taking the union 
of the ranges of N c for all these sets, we find a range for N c which covers the whole range 
for these CKM matrix parameters. We repeat this process for all the form factor models 
mentioned in Eqs. PST) - (1501) and obtain the range of N c for each model as shown in Table 
HTl Taking the union of all the ranges for these models we finally find the maximum possible 
range for N c : 2.74 < N c < 4.77 and 2.81 < N c < 4.92 for q 2 /m 2 b = 0.3 and q 2 /m 2 b = 0.5, 
respectively. 

TABLE II: The range of -/V~ c for all the models and the maximum range of -/V c . 





q 2 /m 2 =0.3 


q 2 /ml=0.5 


Model 1 


(2.74,4.74) 


(2.81,4.90) 


Model 2 


(2.78,4.47) 


(2.84,4.64) 


Model 3 


(2.75,4.77) 


(2.81,4.92) 


Model 4 


(2.76,4.54) 


(2.82,4.72) 


Model 5 


(2.74,4.69) 


(2.81,4.84) 


Model 6 


(2.77,4.50) 


(2.84,4.68) 


maximum range 


(2.74,4.77) 


(2.81,4.92) 



V. SUMMARY AND DISCUSSION 

We have calculated the CP violating asymmetry in the process B° — > p°(uj)p°(uj) — > 
7r + 7T~7r + 7r~ including the effect of p — uj mixing. The advantage of p — uj mixing is that it 
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makes the strong phase difference, 5, between the hadronic matrix elements of the tree and 
penguin operators very large at the u resonance for a fixed N c . We have found that sin<5 
becomes large and reaches the maximum point at the u resonance. Consequently, the CP 
violating asymmetry reaches the maximum value when the invariant masses of the 7r + 7r~ 
pairs in the decay product are in the vicinity of the u resonance. Furthermore, since there 
are two p(u) mesons in the intermediate state, p — uo mixing contributes twice when we work 
to the first order of isospin violation. This leads to an even larger CP violating asymmetry 
than in the case where only single p — uo mixing is involved. This is unique for the process 
B° — > p°(uo)p°(uo) — > 7r + 7r _ 7r + 7r _ . As a result, the largest CP violating asymmetry could be 
more than 90% for some values of N c . This could be observed in the future experiments at 
LHC. Now, we roughly estimate the possibility to observe this CP violating asymmetry. If 
the branching ratio for B° — > p°p° is of order 10~ 6 , then the number of B°B° pairs needed 
for observing the CP violating asymmetry (90%) is roughly BR(B'^—>p p ) a? ~ 

10 6 for la 

signature and 10 for 3a signature 29]. It has been pointed out that at LHC, the number of 
B°B° pairs could be around 4 x 10 7 (for ATLAS and CMS) and 4 x 10 5 (for LHCb) per year 
301 ] . Therefore, it is possible to observe the CP violation for B° — > p°(u)p (u) — > it + -K~it + -n~ 
via the double p — uo mixing mechanism at LHC. 

In the calculation, we need the Wilson coefficients for the tree and penguin operators at 
the scale m^. We work with the renormalization scheme independent Wilson coefficients. We 
have found that apart from the p — uo mixing mechanism, the Wilson coefficients themselves 
could also give observable CP violating asymmetry in some cases. The errors in the CKM 
matrix elements lead to some uncertainties in the CP violating asymmetry. Even bigger 
uncertainties come from the hadronic matrix elements of the tree and penguin operators due 
to the nonperturbative QCD effects. We have worked in the factorization approach, with 
the effective parameter N c being introduced to account for the nonfactorizable effects. We 
have shown that the CP violating asymmetry in this decay process strongly depends on the 
parameter N c . 

In order to determine the range of N c we have compared the theoretical values and the 
experimental data for the branching ratios for B° —>■ p°p° and B° — ► p + p~ . We have found 
that the latest experimental data constrain N c to be in the range (2.74,4.77) for g 2 /m^=0.3 
and (2.81,4.92) for g 2 /m^=0.5, respectively, when we let the CKM matrix elements vary in 
the ranges determined by the current experiments. We have studied the sign of sin5 in the 
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range of N c and found that sind changes its sign at the point N c = 3.68. This also leads 
to the change of the sign of the CP violating asymmetry. Due to the large errors in the 
current experimental data for the branching ratios for B° — > p°p° and B° — > p + p~ we cannot 
constrain N c more accurately at present. If the future experimental data could constrain 
iV c to be less than 3.68 (N c is usually less than 3.68 in other studies), the sign of the CP 
violating asymmetry would remain unchanged in the whole range of N c . Then one could 
remove the mod (ir) ambiguity in the determination of the CP violating phase angle a by 
measuring the CP violating asymmetry in B° — > p°(u)p°(u) — > 7t + 7t~7t + tt~ . 

For the decay process B° — > p°(u)p°(uj), the factorization approach we have used is 
expected to be a good approximation since B meson decays are energetic and since a s (m&) 
and l/m& corrections should be small in the QCD factorization scheme. One may also work 
in the QCD factorization scheme, taking the value of N c to be 3 and including corrections 



of order a s (rrib) as done in Ref. 3jJ. However, the QCD factorization scheme suffers from 
endpoint singularities which are not well controlled. The CP violating asymmetry depends 
on the unknown parameters which are associated with such endpoint singularities. This 



lead to very uncertain CP violating asymmetries in the QCD factorization scheme [311 ] . As 
mentioned before, the uncertainty for the CP violating asymmetry is also very large in the 
factorization approach we have used, i.e. from about -96% to about 56% depending on the 
value of N c and the CKM matrix elements. Furthermore, the CP violating asymmetry may 
strongly depend on the factorization approach adopted 31] . All these issues need further 
and more careful investigations. 



Acknowledgments 

This work was supported in part by National Natural Science Foundation of China 
(Project Number 10675022), the Key Project of Chinese Ministry of Education (Project 
Number 106024) and the Special Grants from Beijing Normal University. 



18 



[1] J. Christenson et al, Phys. Rev. Lett. 13, 138 (1964). 
[2] N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963). 

[3] M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 652 (1973). 

[4] A.B. Carter and A.I. Sanda, Phys. Rev. Lett. 45, 952 (1980); Phys. Rev. D23, 1567 (1981); 

I.I. Bigi and A.I. Sanda, Nucl. Phys. B193, 85 (1981). 
[5] Belle Collaboration, A. Abashian et al, Phys. Rev. Lett. 86, 2509 (2001); BaBar Collabo- 



201 ] and references 



ration, B. Aubert et al., Phys. Rev. Lett. 86, 2515 (2001); see also Ref. 
therein. 

[6] R. Enomoto and M. Tanabashi, Phys. Lett. B386, 413 (1996). 
[7] S. Gardner, H.B. O'Connell, and A.W. Thomas, Phys. Rev. Lett. 80, 1834 (1998). 
[8] X.-H. Guo and A.W. Thomas, Phys. Rev. D58, 096013 (1998); Phys. Rev. D61, 116009 
(2000). 

[9] X.-H. Guo, O. Leitner, and A.W.Thomas, Phys. Rev. D63, 056012 (2001). 
[10] O. Leitner, X.-H. Guo, and A.W. Thomas, Eur. Phys. J. C31, 215 (2003). 
[11] J.D. Bjorken, Nucl. Phys. (Proc. Suppl.) 11, 325 (1989). 
[12] M.J. Dugan and B. Grinstein, Phys. Lett. B255, 583 (1991). 

[13] M. Beneke, G. Buchalla, M. Neubert, and C.T. Sachrajda, Phys. Rev. Lett. 83, 1914 (1999); 

Nucl. Phys. B606, 245 (2001). 
[14] G. Buchalla, A.J. Buras, and M.E. Lautenbacher, Rev. Mod. Phys. 68, 1125 (1996). 
[15] A.J. Buras, M. Jamin, M.E. Lautenbacher, and P.H. Weisz, Nucl. Phys. B370, 69 (1992); 

A.J. Buras, M. Jamin, M.E. Lautenbacher, and P.H. Weisz, Nucl. Phys. B400, 37 (1992); 

A.J. Buras, M. Jamin, M.E. Lautenbacher, and P.H. Weisz, ibid., 75 (1992); M. Ciuchini, E. 

Franco, G. Martinelli, and L. Reina, Nucl. Phys. B415, 403 (1994). 
[16] N.G. Deshpande and X.-G. He, Phys. Rev. Lett. 74, 26 (1995). 

[17] R. Fleischer, Int. J. Mod. Phys. A12, 2459 (1997); Z. Phys. C62, 81 (1994); Z. Phys. C58, 
483 (1993). 

[18] G. Kramer, W. Palmer, and H. Simma, Nucl. Phys. B428, 77 (1994). 

[19] L. Wolfenstein, Phys. Rev. Lett. 51, 1945 (1983); Phys. Rev. Lett. 13, 562 (1964). 

[20] W.-M. Yao et al, J. Phys. G 33, 1 (2006). 



19 



[21] H.B. O'Connell, A.W. Thomas, and A.G. Williams, Nucl. Phys. A623, 559 (1997); K. Malt- 
man, H.B. O'Connell, and A.G. Williams, Phys. Lett. B376, 19 (1996). 

[22] S. Gardner and H.B. O'Connell, Phys. Rev. D57, 2716 (1998). 

[23] Y.-H. Chen, H.-Y. Cheng, B. Tseng, and K.-C. Yang, Phys. Rev. D60, 094014 (1999). 

[24] M. Wirbel, B. Stech, and M. Bauer, Z. Phys. C29, 637 (1985); M. Bauer, B. Stech, and M. 

Wirbel, Z. Phys. C34, 103 (1987). 

[25] X.-H. Guo and T. Huang, Phys. Rev. D43, 2931 (1991). 

[26] D. Melikhov and B. Stech, Phys. Rev. D62, 014006 (2000). 

[27] P. Ball, J. High Energy Phys. 09, 005 (1998). 

[28] P. Ball and V.M. Braun, Phys. Rev. D58, 094016 (1998). 

[29] D. Du, Phys. Rev. D34, 3428 (1986). 

[30] A. Schopper, arXiv: |hep-ex/0605113 

[31] O. Leitner, X.-H. Guo, and A.W. Thomas, J. Phys. G31, 199 (2005). 



20 



Figure Captions 



Fig. 1 The CP violating asymmetry, a, for B° — ► p°{ui)p°{u)) — > 7r + 7r~7r + 7r~. (a): for q 2 /ml = 
0.3, iV c = 2.74 and 4.77 and limiting values of the CKM matrix elements, p and 77: the solid 
line (dotted line) corresponds to N c = 2.74 and maximum (minimum) p and 77; the dashed 
line (dot-dashed line) corresponds to N c = 4.77 and maximum (minimum) p and 77. (b): for 
q 2 jm\ = 0.5, N c = 2.81 and 4.92 and limiting values of p and 77: the solid line (dotted line) 
corresponds to N c = 2.81 and maximum (minimum) p and 77; the dashed line (dot-dashed 
line) corresponds to N c = 4.92 and maximum (minimum) p and 77. 

Fig. 2 Comparison between the CP violating asymmetries in the cases where single and double 
p—oj mixing is involved, respectively, (a): for q 2 /m 2 = 0.3 and N c = 2.74, the solid line (dotted 
line) corresponds to the case with double p — to mixing and maximum (minimum) p and 77; the 
dashed line (dot-dashed line) corresponds to the case with single p — u> mixing and maximum 
(minimum) p and 77. (b): same as (a) but for q 2 /m 2 = 0.5 and iV c = 2.81. 

Fig. 3 sinJ as a function of y/s for B° — > p°(uj)p°(u;) — > tt + tt^tt + tt. (a): for q 2 jm\ = 0.3 and 
N c = 2.74(4.77): the solid line (dashed line) corresponds to the case with double p — io mixing; 
the dotted line (dot-dashed line) corresponds to the case with single p — to mixing, (b): same 
as (a) but for q 2 /m 2 = 0.5 and N c = 2.81(4.92). 

Fig. 4 The ratio of penguin to tree amplitudes, r, as a function of yfs, for limiting values of 
p and m the solid line (dotted line) corresponds to the case of double (single) p — 00 mixing 
with maximum p and 77; the dashed line (dot-dashed line) corresponds to the case of double 
(single) p — to mixing with minimum p and rj. In (a) q 2 /m 2 = 0.3, N c = 2.74 (left) and 4.77 
(right) while in (b) q 2 /m 2 = 0.5, N c = 2.81 (left) and 4.92 (right). 

Fig. 5 The CP violating asymmetry, a, as a function of -/V c and \fs, for B° — > p° (lo) p° (ui) — > 
7r + 7r~7r + 7r for p = p max and 77 = rj max . (a) and (b) correspond to q 2 /m 2 = 0.3 and q 2 /m 2 = 0.5, 
respectively. 

Fig. 6 sinJ as a function of yfs and iV c . (a) and (b) correspond to q 2 /m 2 = 0.3 and q 2 /m 2 = 0.5, 
respectively. 

Fig. 7 Branching ratio for B° — > p°p° for all the models when q 2 /m 2 = 0.3, A = 0.2272, 
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A = 0.818, p = 0.246, and rj = 0.334: the lower (upper) solid line corresponds to Model 1 (2), 
the lower (upper) dotted line corresponds to Model 3 (4) and the lower (upper) dashed line 
corresponds to Model 5 (6). 

Fig. 8 Branching ratio for B° — > p + p~ for all the models when q 2 /mf = 0.3, A = 0.2272, 
A = 0.818, p = 0.246, and -q = 0.334: the lower (upper) solid line corresponds to Model 1 (2), 
the lower (upper) dotted line corresponds to Model 3 (4) and the lower (upper) dashed line 
corresponds to Model 5 (6). 
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Fig. 4(a) 
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Fig. 4(b) 




Fig. 5(a) 
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Fig. 5(b) 




Fig. 6(a) 





c 



Fig. 7 




26 



